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The role of black hole spin and magnetic field threading the 
unstable neutrino disk in Gamma Ray Bursts 
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ABSTRACT 

Aims. We report on the third phase of our study of the neutrino-cooled hyperaccreting torus around a black hole that powers the jet in 
Gamma Ray Bursts. We focus on the influence of the black hole spin on the properties of the torus. 

Methods. The structure of a stationary torus around the Kerr black hole is solved numerically. We take into account the detailed 
treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. 

Results. We find, that in the case of rapidly rotating black holes, the thermal instability discussed in our previous work is enhanced 
and develops for much lower accretion rates. We also find the important role of the energy transfer from the rotating black hole to the 
torus, via the magnetic coupling. 
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1. Introduction 

Observed Gamma Ray bursts (GRBs) form at least two classes 
of events. Long GRBs (Tgg > 2 s) are the the result of a collapsar 
in the core of an exploding massive star, as supported by the de- 
tections of their afterglows. Shorter bursts (Tw < 2 s) are likely 
to be due to the coalescence of a binary system containing two 
compact stars, as observationally supported by the localization 
of their afterglows at the outskirts of their host galaxies (see e.g. 
the review of Zhang 2007). Moreover, it is possible that a third 
class of events exists, namely a fraction of the shortest bursts 
(Tgo < 0.1 s), that can be produced via the evaporation of the 
primordial black holes (Cline et al. 2005). 

Both the two main classes of bursts involve a stage of a newly 
formed black hole surrounded by an accreting torus, responsible 
for the production of a collimated jet. This jet is a source of a 
highly beamed emission observed in the gamma ray band. The 
accreting torus is either fed by the material from the stellar enve- 
lope, which part falls back after the hypernova explosion, or con- 
sumes the material from the disrupted neutron star debris. The 
torus accretion proceeds with highly hyper-Eddington rates, on 
the order of a solar mass per second. Such hyper-accreting tori 
have been already discussed in a number of works (e.g. Popham 
et al. 1999; Di Matteo et al. 2002; Kohri et al. 2005; Chen & 
Beloborodov 2007; Janiuk et al. 2004; 2007; Lei et al. 2009). 

In the present paper we expand on our previous work and 
we further investigate the properties and evolution of such hot 
and dense accretion tori. At the extreme densities and temper- 
atures, determined by the hyper-Eddington accretion rate, the 
torus is cooled mainly by neutrino emission produced primar- 
ily by electron and positron capture on nucleons (f3 reactions). 
The model we develop here was first described in Janiuk et al. 
(2004) for the case of a non-rotating black hole and the nuclear 



matter with a simplified equation of state. We further developed 
this model in Janiuk et al. (2007), considering the more elaborate 
equation of state to account for the microphysics of the accret- 
ing plasma. In the latter work, we solved for the disk structure 
and its time evolution by introducing the new EOS which in- 
cludes photodisintegration of helium, the condition of [5 equi- 
librium, and neutrino opacities. We self-consistently calculated 
the chemical equilibrium in the gas consisting of helium, free 
protons, neutrons, and electron-positron pairs and compute the 
chemical potentials of the species, as well as the electron frac- 
tion throughout the disk. We found an important property of our 
solutions: for very large accretion rates, the torus becomes vis- 
cously and thermally unstable in a narrow range of radii near the 
central black hole. The instability results as the intrinsic property 
of the torus, when the helium is fully photodisintegrated in the 
neutrino opaque plasma. 

In the present work, we further expand our model to account 
for the black hole rotation. The rotating black hole should natu- 
rally be produced in the center of a collapsar or after the merger 
event, when the collapsing material possesses a large amount of 
angular momentum. As a consequence, the black hole spin is a 
plausible mechanism that helps to launch the jets emitted in a 
narrow cone along the rotation axis. As shown in Janiuk et al. 
(2008), especially for the long duration GRBs the BH spin is 
required to sustain the central engine activity and the jet produc- 
tion for a sufficiently long time. Because these long GRBs often 
exhibit very variable temporal structure (e.g. Beloborodov, Stern 
& Svensson 2000), it seems very important to study a possible 
instability mechanism that may account for the variability. We 
check here, for what black hole spin the instability may oper- 
ate near the central black hole, depending on the accretion rate 
and viscosity in the disk. We also check what effects may be im- 
posed on the torus, and its unstable strip, by the additional heat- 
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ing via the magnetic fields, threaded by the rotating black hole. 
The content of the article is as follows. First, we discuss the ba- 
sic assumptions and present the main equations of the model. 
Second, we introduce the changes and corrections to the model 
that are the result of a black hole rotation with an arbitrary spin. 
Third, we introduce the energy extraction from the rotating black 
hole via the magnetic field, as an additional physical process that 
may operate in the gamma ray burst central engine. In Section 
[3]we present the results of our calculations, and in Section|4]we 
discuss the results and conclude. 



2. Model 

The model of a neutrino cooled accreting torus was fully dis- 
cussed in Janiuk et al. (2004; 2007). Here we briefly repeat the 
basic assumptions and the model equations are given in Sec. 
12.11 In Section [2721 we describe the relativistic corrections to the 
model equations, which account for the black hole rotation in 
the Kerr metric. In Section l2~3l we discuss the additional heating 
term in the torus energy balance, that comes from the rapid black 
hole rotation and magnetic dynamo. 

2.1. Basic equations 

The starting model for the torus evolution is calculated by solv- 
ing the energy balance equation at a given accretion rate, M, 
black hole mass, M, and viscosity, a. The steady-state model is 
vertically averaged, with a surface density S = Hp, and height 
H = c s /Qk- Here the sound speed is c s = ^P/p, P is the total 

pressure and Q K = -JGM/r 3 is the Keplerian angular velocity. 
The viscosity is described as in Shakura & Sunyaev (1973), with 
Tnp - -aP and < a < 1 is a model parameter. 

Throughout the calculations we adopt a black hole mass of 
M = 3M Q and several chosen values of a and M. 

For the equation of state, we assume that the torus consists of 
helium, electron-positron pairs, free neutrons and protons. The 
total pressure is due to the free nuclei and pairs, helium, radiation 
and the trapped neutrinos: 



P = P, 



nucl 



^He + ^rad + Py 



(1) 



where f nuc i includes free neutrons, protons, and the electron- 
positron pair gas in beta equilibrium, taking into account the 
partial degeneracy and relativistic temperatures of the species 
(Yuan 2005; Yuan & Heyl 2005). The helium is generally non- 
relativistic and non-degenerate; therefore, its pressure is given 
by the ideal gas EOS. Its number density is calculated using the 
formula of Popham et al. (1999) 



X„ 



295.5p7 3/4 r' 1 /8 exp(-0.8209/r 11 ). 



(2) 



The radiation pressure includes the factor due to the electron- 
positron pairs, and the neutrino pressure is calculated for the 
trapped neutrinos, using the two-stream approximation (Popham 
& Narayan 1995; Di Matteo et al. 2002). We consider both the 
neutrino transparent and opaque regions, as well as the transi- 
tion between the two, and in order to determine the distribution 
function of the partially trapped neutrinos we use a "gray body" 
model (Sawyer 2003), consistent with the two-stream approxi- 
mation. 

The chemical potentials, or equivalently the ratio of free pro- 
tons, x = n p / rib, are determined from the condition of equilib- 
rium between the transition reactions from neutrons to protons 
and from protons to neutrons for a given baryon number density, 



«b, and temperature, T. These reactions are electron and positron 
capture on nucleons, and neutron decay (see Kohri, Narayan and 
Piran 2005; Janiuk et al. 2007). The closing equations for the 
EOS are the conservation of the baryon number, n„+n p = ni,X aac , 
and charge neutrality (Yuan 2005). 

For the neutrino cooling of the torus, we consider the 
electron-positron pair annihilation, bremsstrahlung, plasmon de- 
cay and beta reactions. Each of these emission processes has a 
reverse one, which leads to neutrino absorption, and we calculate 
the absorptive optical depth r a v . for the neutrinos of the three 
flavors. In addition, the free escape of neutrinos from the disc 
is limited by scattering, and we calculate the scattering optical 
depth t s . The neutrino cooling rate is therefore in the neutrino- 
thick torus given by 



^r 4 

3 
4 



T a,yi+Ts 



1 



i=e,(j 



+ -J- + 
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(3) 



Apart from the neutrino emission, the disk is also cooled by 
advection, radiation and photodissociation of helium nuclei. The 
advective cooling is determined approximately as: 



^ dS aPHT 



(4) 



where g a( j v dlnS /dlnr oc (dlnT/dlnr - (F^ - 

l)d In p/d In r) m 1.0. The entropy density S, similarly to the 
total pressure, is the sum of four components: 



S — S nucl + S He + S r ad + S y . 

and we calculate the entropy from the EOS. 
The radiative cooling is 



e rad = 



3P rad c 11 crJ 4 



4r 



AkL 



(5) 



(6) 



where we adopt the Rosseland-mean opacity k = 0.4 + 0.64 x 
10 23 pr~ 3 [cm 2 g _1 ], and the coefficient 11/4 accounts for the 
contribution of electron-positron pairs. 

The cooling rate due to the photodisintegration of a particles 

is: 



Q 



6.28 x 10 28 piov r 



photo 



dX m 



dr 



-H 



(7) 



and X lwc is given by Equation[2] 

Finally, the hydrodynamic equations we solve to calculate 
the disc structure are the standard mass, energy and momen- 
tum conservation. The viscous heating rate is determined by the 
shear, t, v and can be written in terms of a: 



(8) 



The total flux of energy generated at the radius r is deter- 
mined by the global model parameters, i.e. the black hole mass 
and accretion rate: 



3GMM 



(9) 



where f(r) stands for the inner boundary condition and will be 
given below. In order to calculate the initial stationary configu- 
ration, we solve the energy balance: 



Gvisc _ Gadv + 2rad + Qv + 6photo- 



(10) 
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The above equations are appropriate in the Schwarzschild 
black hole case, studied in our previous work (Janiuk et al. 2004; 
2007), and must be modified in the Kerr metric. In Sec. 12.21 we 
include the appropriate correction factors to the disk structure 
equations to introduce the Kerr metric and account for the black 
hole rotation. 



2.1 .1 . Instability in the disk 

The stability analysis of the neutrino-ccoled disk in yS-equlibrium 
was presented presented in Janiuk et al. (2007). The criterion for 
a viscously unstable disk is: 



dM 
5E 



< 



(11) 



and means that, at a given radius, the underdense (overdense) 
region evolves slower (faster), which results in material being 
evacuated from some rings of the disk, while being piled up in 
the others. The unstable region can be conveniently located in 
the steady-state solutions, on the surface density - accretion rate 
(or temperature) plane. On this plane, the branch of thermal equi- 
librium solutions, which has a negative slope, is unstable to the 
surface density perturbations. This branch is also thermally un- 
stable, and we have 



D 



_ 6M*L + 8a ^M^d _ 3 a 2 M 2 



dx 



(16) 



The above coefficients approach unity at large radii and for 
small spin parameters. The inner boundary of the disk is located 
at the last stable circular orbit, depending on the black hole spin 
(Bardeen, Press & Teukolsky 1972): 



GM 



rms = -5- (3 + Z2+ V(3-zi)(3+zi +2z 2 )) 



(17) 



where zi = l+(l-a 2 ) 1/3 ((l+a) 1/3 + (l-a) 1/3 ),z 2 = (3a 2 +z 2 ) 1/2 , 
and the upper sign in EqfTT] is for a retrograde orbit while the 
bottom sign is for direct orbit around a black hole. (In this work 
we always consider the direct orbits.) 

The increase of gravity in the close vicinity of the rotating 
black hole in consequence leads to a smaller disk height, which 
will be given by the hydrostatic equilibrium: 



(18) 



The angular velocity in the disk around a spinning black hole 
is given by: 




31og0 + 3 logg- 



ed log T d log T 



(12) 



Q = 



and therefore any small increase (decrease) in temperature leads 
to more (less) heating rate with respect to the cooling. 

We found, that for a non-rotating black hole, the neutrino- 
cooled disk was generally stable for all radii if M < 10M o s , 
while at larger accretion rates the unstable region appears at 
the innermost radii. The instability grows, because the energy 
balance changes because of the photodisintegration of helium, 
while the neutrinos are trapped in the opaque plasma. In the un- 
stable region, the electrons and protons become non-degenerate 
and the electron fraction rises inwards in the disk. The changes in 
the electron fraction influence the total pressure and the energy 
energy dissipation rate, so that the heating and cooling balance 
is affected and the local accretion rate increases. We found that 
for the highest accretion rates we used (12 M s _1 ) and the non- 
rotating black hole, the radial extension of the unstable part of 
the disk was up to ~ 8 Schwarzschild radii. 

2.2. Black hole rotation 

The structure of the relativistic accretion disk is calculated using 
the correction factors derived in Riffert & Herold (1995) for a 
Kerr black hole (see e.g. Reynoso, Romero & Sampayo 2006). 
The black hole spin is parameterized by a dimensionless specific 
angular momentum (Bardeen 1970), a = cJ/GM 2 , and the fol- 
lowing functions of a, M, and radial coordinate r are introduced: 



2GM l GMa \ 2 

c 2 r \ c 2 r I 



5=1- 



3GM „ (GM\ V2 
+ 2a — 



(GM\ 3/2 I GMa \ 2 



(13) 
(14) 
(15) 



GM 



The viscous shear will be modified as: 
A 

Tnp = ~aP- 



BC 



(19) 



(20) 



Consequently, the heating and cooling terms in the energy bal- 
ance equation will be affected by the above correction factors. 
Also, the equation for the energy flux generated in the disk, Eq. 
©, will now be specified in the Kerr metric. Therefore, instead 
of the Newtonian boundary condition, we use the function: 



B 



(21) 



This factor is equal to zero at the r ms radius and approaches 
unity at large radii. This asymptotic behaviour of f(r) is the same 
as for the boundary condition derived in Chen & Beloborodov 
(2007), who used more complex formalism in the Kerr metric. 
The torque imposed on the disk due to magnetic coupling with 
the black hole (see below) does not change the no-torque inner 
boundary condition either. 

2.3. Magnetic heating 

The spinning black hole may transfer its rotation energy to the 
jet via the Blandford-Znajek (BZ) process (Blandford & Znajek 
1977), as well as to the accretion disk, via the magnetic coupling 
(Li & Paczyriski 2000). In the latter, the closed lines of magnetic 
field connect the black hole horizon with the disk, and if the 
BH rotates faster than the surrounding disk, the torque is exerted 
and the energy and angular momentum are extracted from the 
BH. On the other hand, the direction of energy transfer may be 
reversed, if the disk angular velocity is larger than that of the 
black hole. 

The details of the process depend on the structure of mag- 
netic field. Here we follow the assumptions of Wang et al. 
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(2003), that the magnetic field is constant at the horizon and 
varies as a power law with the disk radius. We also use the map- 
ping relation between the angular coordinate, 0, on the BH hori- 
zon and the disk radius, r, derived by these authors. Due to this 
relation, the closed magnetic filed lines are attached not only 
to the inner edge of the disk, but connect the horizon region of 
smaller 6 with outer disk radii, up to some maximum radius r out . 
Therefore, the disk inner ring, between r ms and r out , will be af- 
fected by the energy transfer to and from the rotating black hole, 
depending on its spin and magnetic field configuration. Because 
the inner regions of the disk are also subject to the thermal in- 
stability, we intend to check what effect on this instability may 
arise due to the energy transfer between the black hole and the 
disk. 

The mapping relation is defined by the following change of 
variables (Wang et al. 2003): 



sin 6d6 - -G(a, ri)d^ 



and 



cos 6 ■ 



where 



r 



G(a,i;,n)dt; 



(22) 



(23) 



G(a,i;,n) 



2 V(l + « 2 * m 4 s + 2^X1 - 2 X £f l + a^xlt 1 ) 



Here the dimensionless radius is £ = r/r ms , x — V r / r g an d 
r g = GM/c 2 . The mapping function G(a,£,n) depends on the 
black hole spin, which is a parameter of our model, and on the 
strength of magnetic field at a given radius. For the magnetic 
field, is assumed that it varies as (Blandford 1977): 



B z oc if" 



(25) 



and the index « is a model parameter. 

The magnetic coupling between the rotating black hole and 
the disk exerts a torque, equals to: 



r-K/2 

Tuc=T 4a(\+ Vl -a 1 ) ^ 



(1 -frsvcfGdO 
2-(l- Vl - a 2 ) sin 2 6 



, (26) 



where B = Q/Q# is the ratio between the angular velocity in the 
disk, given by Eq.QjO and the velocity at the black hole horizon, 
given by 



Oh =al(lr g {\ + Vl - a 2 )) . 



2 -2 

cm" s - 



In general, the proportionality coefficient in the above equation 
can be on the order of 1.0, however, some departures from the 
equipartition of magnetic energy are possible. In case of the a- 
disks, the ratio of B 2 1 /87rf max = B m . dg cannot be too high, and will 
rather be on the order of a (Moderski, private communication). 

Finally, the energy transfer from the rotating black hole to 
the accretion disk results in the additional heating term: 



(30) 



(27) 



The integration over coordinate 9 from 8 m i n to # max is equivalent 
to the integration over disk radius from £ max to £ m i n . Here we set 
£min = £ms = 1 at 6 — n/2 and for simplicity we assume that 
no open magnetic field lines take the black hole rotation power 
to infinity, i.e. the BZ effect is vanishing and £ max has its largest 
value for 6 m [„ = 0. 

The normalization of the torque is (Wang et al. 2003): 



(28) 



and the energy density of the magnetic field is proportional to 
the total pressure in the disk (Moderski, Sikora & Lasota 1998): 



B\ oc 87rP max . 



(29) 



3. Results 

3. 1 . Disk instability for a rotating black hole 

We found that the thermal instability of the neutrino cooled disk, 
previously studied for the case of Schwarzschild black holes in 
Janiuk et al. (2007), is also present in the disks of Kerr black 
holes. The instability arises in the inner parts of the disk, where 
the helium is fully photodisintegrated, the nuclei pressure is lo- 
cally lower, while the neutrino opacity is large and the dominant 
pressure terms are due to the neutrinos and radiation. In our pre- 
vious work, we found for a non-rotating black hole that the radial 
extension of the unstable strip was up to about 8 Schwarzschild 
radii. The accretion rate for which the instability appeared, had 
to be extremely large, on the order of 12 solar masses per sec- 
, 9 ond. The instability was manifested in both stationary models, 
as a phase transition to a distinct branch of distribution for vari- 
ous quantities (temperature, density, electron fraction, chemical 
potentials of species, etc.), as well as in the time dependent evo- 
lution of the disk (variability in time). 

In this work, we calculated a grid of stationary models in 
the Kerr metric, to check for the location of the unstable strip 
depending on the black hole spin. We used moderate accretion 
rates and several viscosity parameters. In Figure [TJ we show the 
radial extension of the unstable strip, as a function of spin pa- 
rameter a, for three values of accretion rate: M = 0.5, 1.0 and 
5.0 M Q s _1 . We find that for the extreme Kerr black hole, the 
instability appears already at very low accretion rates, 0.5 solar 
masses per second. The dependence on a is rather weak. In gen- 
eral, the larger is viscosity, the larger BH spin is needed for the 
thermal instability to appear at a given accretion rate. Also, for 
larger alpha, the instability strip is narrower. As shown in the up- 
per panel of the Figure, only for a > 0.95 at viscosity a = 0.1, 
and for a > 0.98 at a = 0.3 the unstable solutions are found 
at this low accretion rate. For lower accretion rates, we checked 
M — 0.1 M s" 1 , and no unstable solutions were found down to 
a = 0.03 and at a = 0.998. 

In general, the predicted trend is conserved: the higher the 
accretion rate and the larger black hole spin, the larger is the 
extension of the unstable strip. The trend with viscosity is the 
following: for larger viscosity parameters the extension of the 
unstable zone shrinks. 

In the Figure [2] we show the dominant cooling terms as a 
function of radius. The terms are due to the neutrino emission, 
Q v , and helium photodisintegration, Qphoto- The other cooling 
terms, due to radiation and advection of energy, are incorporated 
in our model, but they are a few orders of magnitude smaller 
than neutrino cooling. Therefore we do not show these terms on 
the plot. The 2 p h to term at large radii can locally have a nega- 
tive value, which means that the helium nuclei are at this radius 
created. 

The unstable strip in the accretion disk can be localized be- 
tween the two spikes in the Qphoto distribution. The first spike is 
very close to the inner edge of the accretion disk, with a small 
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Fig. 1. The parameter space for which the unstable solution is 
found, in the plane radius - spin parameter. The bottom panel 
shows the results for accretion rate M — 5.0 M Q s _1 , the middle 
panel is for M — 1.0 M s and the top panel is for M = 0.5 
M s . The two values of viscosity were used: a = 0.1 (shaded 
area) and a = 0.3 (cross-shaded area). 



15 







ii 


1 

m=i.o : 

a = 0.95, 0.98 - 




/ 

1 ^ 

1 ,' 






L- 






— 1 1 1 




1 (- 


i 1 1 — 

VQ, 


\ 

M = 0.5 : 
a = 0.95, 0.98 1 






Qphoto 




^ V 



0.5 1 
log r/r g 



Fig. 2. The main cooling terms in the disk as a function of radius: 
neutrino cooling (thick blue lines) and helium photodisintegra- 
tion (thin red lines). The bottom panel shows the results for ac- 
cretion rate M — 0.5 M Q s , and the upper panel is for M — 1 .0 
M s _1 . The two values of black hole spin parameter are shown: 
a = 0.98 (solid lines) and a - 0.95 (dashed lines). The viscosity 
parameter is a = 0. 1 . 



offset because of the inner boundary condition (cf. Eq. [2TT>. Due 
to this condition, no torque is induced on the inner edge. This 
results in zero viscous dissipation rate and the plasma temper- 
ature drops close to the inner edge (in practice, for numerical 
reasons we set our inner radius of the computational domain 
slightly above inner edge, to avoid singularities). Note, that the 
torque induced by the magnetic coupling between the disk and 
the black hole is transported only outwards with radius, while 
at the inner edge of the disk the torque is still zero. We also 
note, that the density drop at the inner edge is compensated by a 
large radial velocity, because we take into account the advective 
cooling. Therefore, despite the small disk thickness, the accre- 
tion rate is kept constant. Above the inner edge of the torus, the 
temperature reaches 10 10 - 10 11 K, while the density exceeds 
10 9 g cm -3 , and the helium nuclei can appear. The plasma con- 
sists of free neutrons, protons and electron-positron pairs, and if 
the photon flux is sufficient, the helium nuclei are photodisinte- 
grated. Depending on the model parameters, the outer radius of 
this strip can reach up to 12 gravitational radii (Fig.[TJ. Outwards 
of this radius, the helium is again synthesized and its photodis- 
integration contributes to the energy balance in the outer disk 
(with some small fluctuations in Qphoto up to r ou t, not visible in 
the Figure). 

The neutrino cooling term is always a dominant one in all the 
models. The larger accretion rate and black hole spin, the hotter 
is the disk and hence the neutrino cooling is larger. For the ro- 
tating black hole, the inner disk radius shifts closer (cf. Eq. [T7J, 
and also the maximum in the heating and cooling distribution is 
shifted inwards. 

We note, that our model works well in the neutrino-cooling 
dominated region, which is always inside the radius r max . This 
'neutrino ignition radius' strongly depends on the accretion rate 
and weakly depends on the black hole spin. Outside this region, 
the flow is dominated by advection and neutrino emission is in- 
efficient in cooling the disk (see also Chen & Beloborodov 2007; 
note however that we do not consider the gravitational instabil- 
ity which according to these authors is present at large radii and 




1.6 - 



-1 -0.5 0.5 

log M [Me s-'] 

Fig. 3. The parameter space for which the disk is dominated by 
the neutrino cooling (NDAF). We plot the maximum radius of 
NDAF (in the units of r g = GM/c 2 ) as a function of accretion 
rate, for 2 chosen values of black hole spin: a=0 (solid line) and 
a=0.9 (dashed line). 



large accretion rates, nor we study the case of the smallest accre- 
tion rates, for which the neutrino cooling is inefficient through- 
out the disk at all radii). In Figure [3] we indicate the parameter 
space for which the neutrino cooling is dominant. The instability 
we discuss here, as well as the magnetically coupled region we 
discuss below, are both located well inside the neutrino cooled 
region. 
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3.2. Magnetic coupling between the disk and rotating black 
hole 

The mapping relation between the polar coordinate on the black 
hole horizon and the disk radius implies that there exists a maxi- 
mum radius, defining the region of influence of the rotating black 
hole. Below this radius, at £ < £ out , the energy may be transferred 
to and/or from the black hole to the disk. Above this radius, there 
are no closed magnetic field lines, and the disk is no longer sub- 
ject to the heating by the black hole energy extraction. 

In the Figure |4] we show the dependence of the radial exten- 
sion of the magnetically coupled region on the black hole spin 
and magnetic field distribution. The latter is parameterized by 
the index n in Eq.[25] and the Figure shows several examples for 
n=1.5, 2.0, 3.0 and 6.0. As was found by Wang et al. (2003), the 
influence of the black hole rotation is concentrated in the closest 
vicinity of the inner disk edge and the radius £ out = r out /r ms is 
very small for small values of index n. For n < 2.0 the largest 
£out is about 2.2, at the black hole spin a = 0.98 and then drops. 
The result for larger n are quantitatively the same as in Wang et 
al. (2003) and they differ only quantitatively. For n — 3 the maxi- 
mum is £ out > 10.0 at a — 0.998, while those authors had a value 
about £ out « 4.0, which in turn we obtained for n = 2.7. The 
reason is probably the sensitivity of the results to the numerical 
method used to compute the integrals in Eq.|26] For n > 3.0 and 
large a, the size of the region £ out goes to infinity, and the larger 
n, the smaller is the critical spin parameter at which that happens 
(see Wang et al. 2002). 

Note that actually in Fig. [4] we plot the results in the units 
of the gravitational radius, r g , instead of r ms radius, i.e. the de- 
pendence on the black hole spin is accounted for also in the r ms 
value. Therefore we can directly compare the extension of the 
magnetically coupled region with the extension of the thermal 
instability in the disk, presented in Fig.Q] In these units, the size 
of the magnetically coupled region decreases with a, because 
also the r ms radius is smaller (inner edge is closer to the black 
hole) if the spin is large. For the most interesting cases of the 
large spin, the outer radius of magnetically coupled region is 
usually well inside the thermal instability strip, provided the ac- 
cretion rate is sufficient for the thermal instability to occur, and 
n is rather small. However if n > 3, the magnetically coupled re- 
gion may extend to larger radii than the thermally unstable strip, 
especially for moderately small accretion rates. 

However the size of the magnetically threaded strip may be 
small, the energy input given to the disk by a rotating black hole 
is quite large. In the Figure [5] we show the maximum magnetic 
torque induced in the disk by the black hole rotation, depend- 
ing on its spin, for several values of index n. The value of the 
torque increases with radius and saturates at £ out at its maxi- 
mum value. This value is correlated to both index n and spin, 
decreasing only for the extreme Kerr black holes. The value of 
maximum Tmc = 0, obtained for the most slowly rotating black 
holes, means that in fact for these parameters the torque is neg- 
ative throughout the magnetically coupled region. This means 
that the rotation energy is transferred only from the disk to the 
black hole. On the other hand, for the fast rotating black holes, 
and large n, the value of 7mc may be large and positive at the 
outer parts of the magnetically coupled region, while negative in 
the innermost radii. This is because if r ms is sufficiently small, 
the innermost disk parts rotate faster than the black hole, and the 
energy is extracted from the inner disk to spin up the black hole. 

Now, we study how the magnetic coupling between the disk 
and the rotating black hole will influence the thermal-viscous in- 
stability, found by Janiuk et al. (2007). This thermal instability 




Fig. 4. The radial extension of the region threaded by magnetic 
coupling between the disk and rotating black hole, as a function 
of its spin. The solid line marks the power law index of magnetic 
field n — 1 .5, the dotted line is for n = 2, long-dashed line is for 
n = 2.5, the ditted-dashed line is for n = 2.7 and the short dashed 
line is for n — 3. 
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Fig. 5. The maximum magnetic torque induced in the disk, as a 
function of the black hole spin. The solid line marks the power 
law index of magnetic field n = 1.5, the dotted line is for n = 2, 
the long dashed line is for n = 2.5, the dotted-dashed line is for 
n = 2.7 and the short dashed line is for n = 3. 



was studied in case of the non-rotating black hole and found to 
require very high accretion rates. In the present work we con- 
firmed the instability occurs also for the rotating black holes and 
needs smaller accretion rates (Section [3~TT i. To further investigate 
the role of the black hole rotation, we calculated a grid of models 
parameterized by the accretion rate M, black hole spin a, as well 
as the magnetic field power law index n and its normalization 
/Wg, which describe the strength of the magnetic coupling. 

In Figure [6] we show the density profiles for two values of 
accretion rate, M = 0.5 and 1.0 M s _1 , and several spin param- 
eters in the range from a = 0.9 to 0.998. We concentrate on the 
innermost region of the accretion disk, where both the thermal- 
viscous instability and the magnetic coupling overlap. Therefore 
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the Figures show the region up to the radius r = 15r ms . As dis- 
cussed in Janiuk et al (2007), the thermal-viscous instability re- 
sults in a local density decrease accompanied by the tempera- 
ture rise in the inner strip of the disk. Outwards of that strip, a 
stable solution is found, for which the density and temperature 
smoothly decrease with radius. 

Now, we compare our previous unstable solutions, calcu- 
lated under the assumption of no magnetic torque (i.e. f3 mdg = 0, 
shown by the thinner lines), with the new solutions with the mag- 
netic torque included (i.e. /3 mag = 0.5 or 0. 1 , shown by the thicker 
lines). We note that the density dip is deeper, whenever the en- 
ergy is transferred from the black hole to the disk. Only very 
close to the inner edge of the disk, for large BH spins and index 
n ~ 3, we have the black hole extracting energy from the disk. 
When this is coupled to the thermal instability, the density at 
the innermost radii increases, and the instability strip effectively 
shifts towards larger radii. For the outer edge of the instability 
strip, we can identify it with the local maximum of the density. 
Therefore, this point is shifted towards a smaller radius and in 
turn the whole unstable strip is narrower. However, if both a and 
n are sufficiently large, the magnetically coupled region extends 
well outside the instability strip (cf. Fig.|6j bottom panel). In this 
case, the density profile is still affected, and the density is sys- 
tematically lower than in the respective model without the mag- 
netic torque. However, now the density also decreases smoothly 
with radius, and we do not find this region unstable. 

For smaller index, n = 2, the magnetically coupled region 
does not cover the whole instability strip. This is shown in the 
top panel of the Figure [6] Therefore we have the density de- 
crease only in some part of that strip, while outwards the solu- 
tion matches the Qmc~ - case. The effect is in these models 
smaller than in case of n — 3, because we now use the y8 mag =0.1 
instead of f} nYd% = 0.5, so that the magnetic heating does not ex- 
ceed the total emitted flux F toi . For n — 2 and a < 0.98, the 
magnetic torque is everywhere positive, i.e. there is no region 
where the energy is extracted from the disk towards the black 
hole, and there is no density rise on the inner edge. 

The local density dip is produced by the magnetic term also 
in these models, where there is no thermal instability. This is 
shown in Fig. [6] for M — 0.5, a = 0.95 (bottom panel) and for 
M = 1.0, a — 0.9 (upper panel). The radial extension, depth and 
localization of this dip are determined by the parameters n, a and 

Anag- 

In the Figure|7]we plot the temperature profiles for the same 
models as were shown in Fig. [6] As we mentioned above, the 
local density dip is accompanied by the temperature excess in 
the unstable region. However, when that region is magnetically 
threaded by a rotating black hole, the excess is smaller and ex- 
tends to a smaller radius than for Qmc = 0. In contrast, the stable 
disk outside the instability strip, is heated up and the temperature 
there increases, because the extra heating term is added. On the 
other hand, the temperature decreases close to the inner edge in 
the models with a negative torque, i.e. when the energy is trans- 
ferred from the disk to the black hole. The stable disk affected 
by the magnetic coupling is therefore cooled at the radii where 
Gmc < and heated at the region of Qmc > 0. 

Such a behaviour of the temperature profile, which shows 
the net decrease of the temperature when the magnetic heating 
is added in the unstable region, is caused by the nature of the 
thermal/viscous instability. As discussed in detail in Janiuk et 
al. (2007), the instability is a result of a phase transition, which 
occurs under certain conditions in the nuclear matter. As a result, 
the extra energy provided by the rotation of the black hole via 
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Fig. 6. The density profiles for the magnetically heated/cooled 
disks (thick lines), for several values of the black hole spin: a = 
0.998 (dot-dashed lines), a = 0.98 (solid lines), a = 0.95 (short 
dashed lines) and a = 0.9 (long dashed lines). The thinner lines 
show the case when the magnetic torque was neglected. The pa- 
rameters were M = 1.0M o s~' (top panel) and M = OSMqS' 1 
(bottom panel). The magnetic torque parameters were n — 2, 
ySmag = 0.1 (top panel), and n = 3, y6 mag = 0.5 (bottom panel). 
The viscosity parameter is a — 0.1. 
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Fig. 7. The temperature profiles for the magnetically 
heated/cooled disks (thick lines), for several values of the 
black hole spin: a = 0.998 (dot-dashed lines), a - 0.98 (solid 
lines), a = 0.95 (short dashed lines) and a - 0.9 (long dashed 
lines). The thinner lines show the case when the magnetic 
torque was neglected. The parameters were M = 1.0M o s _I (top 
panel) and M = 0.5M o s~' (bottom panel). The magnetic torque 
parameters were n — 2, yS mag =0.1 (top panel), and n - 3, 
Anag = 0.5 (bottom panel). The viscosity parameter is a - 0.1. 



the magnetic coupling is consumed by this phase transition, and 
cannot heat the disk. 

We conclude, that the magnetic coupling between the spin- 
ning black hole and the disk does influence the thermal-viscous 
instability. It can result in shifting the unstable region towards 
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larger inner radius and shrinking it in size. The strength of the 
thermal-viscous instability will be modified under the influence 
of the magnetic coupling, because the density dip in the unsta- 
ble region will be deeper. However, the temperature excess in 
the unstable region is smaller, if the disk magnetically coupled 
to the black hole. To quantify the net effect and calculate the re- 
sulting luminosity fluctuations we would need to make the time- 
dependent modeling, which is beyond the scope of the present 
paper. We intend to study this in more detail in the future work. 

We also note that there is some small impact of magnetic 
coupling on the stable disk solutions. This is an extra dip in den- 
sity accompanied with the temperature excess, which arises in 
the previously stable region. However this effect seems to be 
much less pronounced than the thermal-viscous instability. This 
may potentially lead to the development of a separate type of a 
magnetic instability. 

4. Discussion 

The importance of black hole spin in the accretion disk systems 
and the formation of jets, also in the presence of the magnetic 
fields, has been studied by many authors. It has been discussed 
e.g. in the review by Blandford (1999). In the standard accretion 
disk theory (Shakura & Sunyaev 1973) it is assumed that there 
is no coupling between the black hole and the disk. However, if 
the magnetic field is considered, such a coupling may be present. 
Contrary to the Blandford-Znajek process (in which the mag- 
netic filed lines threading the black hole may close on a very 
remote load, far away from the black hole), the magnetic field 
lines coupled to the accretion disk are closed in the vicinity of 
the inner radius. In this process, the energy and angular momen- 
tum are extracted from the black hole and transferred to the disk, 
if the hole rotates faster than the disk (MacDonald & Thorne 
1982; Li 2000). Such coupling may have much higher efficiency 
in extracting the black hole rotation energy than the Blandford- 
Znajek mechanism, as was first quantitatively estimated by Li 
(2002). 

Let us also remind here after Li (2002), that in this model the 
torque produced by the magnetic coupling with the black hole 
propagates only outwards, and no torque is imposed at the inner 
boundary of the disk. The non-zero torque at the inner boundary 
is possible if the disk is magnetically connected with the plasma 
in the transition region between the disk and the black hole (e.g. 
Krolik 1999), which we do not consider here. Another impor- 
tant assumption behind the present model is that about the quasi 
steady-state of the disk and magnetic field configuration. Under 
this assumption, the magnetic field lines are not tangled by the 
rotation and MHD instabilities of the Balbus & Hawley (1991) 
are not taken onto account. In our present work we hold this as- 
sumption, after Li (2002) and Wang et al. (2002; 2003a,b). The 
evolution of the coronal magnetic field threading the differen- 
tially rotating disk was studied however in the literature, e.g. by 
Lovelace et al. (2002) in the context of Poynting flux dominated 
jets in quasars, and in Uzdensky (2004; 2005) by numerically 
solving the Grad-Shafranov equation in the Schwarzschild and 
Kerr metric. 

In the context of Gamma Ray Bursts, the black hole spin in- 
teraction with the magnetized torus was studied by van Putten 
et al. (2004) These authors take into account the dissipation of 
energy in the torus not only via the thermal and neutrino emis- 
sion, but also through the gravitational waves and winds. They 
estimate the lifetime of the central engine powered by rapidly 
spinning black hole, whose spin-down is governed by the sur- 
rounding magnetic field coupled to the inner torus. In the present 



model, we neglect the gravitational waves emission and we con- 
centrate on the structure of the neutrino-cooled disk. We restrict 
ourselves to the steady-state solution, however we note that the 
lifetime of the long gamma ray burst central engine is defined 
by the period of the rapid spin of the black hole, which may de- 
crease with time during the engine activity (see also Janiuk et al. 
2008). 

In this article, we study the model of a stationary neutrino 
cooled disk around a rotating Kerr black hole. We confirmed 
the presence of the thermal-viscous instability in the region of 
the large neutrino opacity and efficient helium photodisintegra- 
tion. The instability was first studied in Janiuk et al. (2007) in 
the time-dependent model for the Schwarzschild black hole. It 
was found to produce dramatic changes in the density and tem- 
perature profiles, as well as a rapid time variability of the disk 
emission. Here, we quantitatively described the role of the black 
hole spin and viscosity in the disk. We found that in the Kerr 
black hole disks the large spin parameter results in a larger size 
of an unstable strip, which appears even for a moderate accretion 
rate. The critical accretion rate is anticorrelated with the BH spin 
required for the onset of the instability. In the extreme BH case, 
the instability close to the inner edge may be effective even for 
M = 0.5M Q s _1 . On the other hand, the more viscous is the disk, 
the smaller in size is the unstable strip and the weaker should be 
the instability. 

The time dependent calculations were not performed in the 
present work, because they would require to compute time evo- 
lution of the BH spin and a moving radial grid, which is intended 
to be studied in the future work. However, we studied the influ- 
ence of the black hole rotation on the disk structure considering 
our stationary model. This was done by means of the magnetic 
transfer of energy from the rotating black hole to the disk, and 
vice versa. 

The first important finding in this work is that thermal in- 
stability occurs for quite low accretion rates, on the order of a 
fraction of solar mass per second, provided the black hole ro- 
tates very fast. This can be relevant for the origin of the variable 
energy input to the GRB jets not only in case of the short GRBs, 
but also for long ones. The latter are presumably powered by the 
accretion of the fallback material from the massive star envelope 
onto the newly born black hole after the hypernova explosion 
(e.g. Mac Fadyen & Woosley 1999). In these models, the accre- 
tion rate is found to be substantially lower than in the merger 
of two neutron stars, and therefore the event can last for much 
longer time to power the GRB for several hundreds of seconds. 
On the other hand, the massive star that is a progenitor of a hy- 
pernova is a rotating star and therefore a natural consequence 
of the explosion should be the rapidly rotating black hole in the 
center. 

As was recently shown by Janiuk, Moderski & Proga (2008), 
the large black hole spin is a key ingredient for the occurrence of 
the longest duration GRBs, and is required for the GRB central 
engine to operate for a time of the order of hundreds of seconds. 
The present study leads us therefore to a general conclusion that 
the longest GRBs, that are powered first by the neutrino anni- 
hilation and later by the black hole spin and require a > 0.9, 
should be very variable at the initial phase of the prompt emis- 
sion. Later on, when the accretion rate drops, the variability may 
be suppressed and what is observed is a smooth tail in the gamma 
ray lightcurve. On the other hand, the period of the highly vari- 
able GRB emission may be somewhat extended if the black hole 
is being spun up by the accretion. The issue discussed in the lit- 
erature was whether in such a case the spin equilibrium value 
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can be achieved at a — 0.998 (Thorne 1974), or a less value 
(Gammie, Shapiro & McKinney 2004). 

In our study, we also considered the heating of the inner- 
most regions of the torus due to the magnetic coupling with the 
rotating black hole. We found that this process can have some in- 
fluence on the thermal instability, because the size and position 
of the unstable region are somewhat changed. Also, the profiles 
of density and temperature in the unstable region are modified. 
The magnetic coupling enhances the change in density profile, 
however it weakens the change in temperature profile. The en- 
ergy transfer to the disk from the rotating black hole does not 
suppress the thermal instability. This finding also supports our 
conclusion that the variable and long duration GRBs are pow- 
ered by the rapidly spinning black holes. 
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